Generalized Ulam-Hyers-Rassias stability and novel sustainable techniques for dynamical analysis of global warming impact on ecosystem

Marine structure changes as a result of climate change, with potential biological implications for human societies and marine ecosystems. These changes include changes in temperatures, flow, discrimination, nutritional inputs, oxygen availability, and acidification of the ocean. In this study, a fractional-order model is constructed using the Caputo fractional operator, which singular and nol-local kernel. A model examines the effects of accelerating global warming on aquatic ecosystems while taking into account variables that change over time, such as the environment and organisms. The positively invariant area also demonstrates positive, bounded solutions of the model treated. The equilibrium states for the occurrence and extinction of fish populations are derived for a feasible solution of the system. We also used fixed-point theorems to analyze the existence and uniqueness of the model. The generalized Ulam-Hyers-Rassias function is used to analyze the stability of the system. To study the impact of the fractional operator through computational simulations, results are generated employing a two-step Lagrange polynomial in the generalized version for the power law kernel and also compared the results with an exponential law and Mittag Leffler kernel. We also produce graphs of the model at various fractional derivative orders to illustrate the important influence that the fractional order has on the different classes of the model with the memory effects of the fractional operator. To help with the oversight of fisheries, this research builds mathematical connections between the natural world and aquatic ecosystems.

through photosynthesis.This biological mechanism is key to comprehending the prospects of environmental change as well as is an important part of the natural carbon cycle on the planet's surface 3 .Temperature variations throughout the year, water column interaction, supply, and consumption all have an impact on how they behave.These variables can be modified by climate, changing the taxonomic composition, seasonal dynamics, and phytoplankton structure.Climate affects phytoplankton in two ways: directly by physiological processes and indirectly by the supply of nutrients, stratification of the water column, and heterotrophic feeding.These changes affect numerous processes 4 .The majority of carbon dioxide and oxygen are produced and absorbed by marine ecosystems.High levels of greenhouse gases cause the oceans to warm quickly and are contributing to global warming.Due to the disruption of the bicarbonate buffer, which keeps ocean acidity between 7.5 and 8.4, this harm has an impact on planktonic species and fisheries resources.Rapid global warming is raising atmospheric carbon dioxide levels, creating an ocean climate that is unfavorable and resulting in decreased plankton and fish populations in marine ecosystems.If this keeps up, by the end of this century a sizable section of marine ecosystems might be eliminated or degraded 5 .In the past 60 years, the Indian Ocean's marine fish and phytoplankton populations have declined by 20% due to global warming.If the problem remains, fish and plankton populations may drop to 50-90% of their current levels, creating an ecological desert and lowering production 6 .By the end of the twenty-first century, the tropical Pacific might experience a temperature increase of more than 3 • C, which would threaten between 50 and 80 percent of marine species, especially plankton 7 .
The movement of people from rural regions is speeding up urbanization trends, causing population density to rise and altering the quality of life globally.This results in environmental contamination, which harms the land, water, and air while ruining the natural world and changing the climate in urban areas 8 .Global environmental changes are a result of human activity, which started with the Industrial Revolution in the 1750s.Before industrialization, carbon dioxide emissions were minimal, but the end of the Ice Age and the Industrial Revolution resulted in a considerable rise in greenhouse gas concentrations 9 .Carbon dioxide storage and capture are potential ways to stop ocean acidification and global warming.On the other hand, leakage from storage structures can hasten acidification, which might have an impact on environmental bacteria.To comprehend the impacts of carbon dioxide, a review of how various bacteria are affected by carbon dioxide is required 10 .The environment is a battleground for living things, and after the beginning of the industrial era, emissions of greenhouse gases have considerably increased.Since the last 0.8 million years, the average atmospheric carbon dioxide concentration has increased from 280.01 to above 380 ppmv, surpassing that record 11 .The most distant continent, Antarctica, is subject to adverse effects from human activities, such as excessive harvesting, environmental impact, and the invasion of alien species.If these problems aren't addressed at the same time, the marine ecosystems in Antarctica will deteriorate and become more similar to other marine ecosystems in terms of substance, structure, and mechanisms 12 .Warming elevations and severe weather events lead to epidemics, even if vector-borne infectious illnesses may not have the most negative effects on health.These factors, such as frequent rainstorms, rodent infestations, and tainted water, render mosquito-borne parasite and viral infections climate-sensitive 13 .
The consequences of global changes in climate on coastal biodiversity and marine ecosystems are extensively studied in the literature 14,15 , with multiple publications 16 statistically detailing possible effects on marine ecosystems and coastline aquatic organisms.Mathematical modeling and esoteric mathematics are included in the field of mathematics.With the use of mathematical ideas and hypotheses, it is simple to assess the progress of the task, the procedure, the forecasts, and the results.As a result, scientists rely heavily on mathematics today 17 .In several publications, mathematical modeling is used to examine how global warming would affect marine ecosystems 18 offered a stochastic mathematical model, while 1,19 proposed a deterministic mathematical model to explain the effects of rising temperatures on aquatic ecosystems.Some papers' writers used statistical-based evaluation to carry out their research 16 , while others used literature-based assessment to demonstrate their findings 20 .
The capacity of fractional calculus to analyze genetic changes and their effects on the dynamic operations of physical structures has drawn attention.While the fractional order concerns include integration and contour differentiation, these non-local aspects help us grasp practical issues like memory characterization and genetic qualities 21 .Numerous fractional operators are employed in the mathematical modelling of problems that arise in the actual world [22][23][24][25] .Sekerci and Ozarslan 26 investigated the consequences of predation on the oxygenplankton system by analyzing a fractional model of oxygen, phytoplankton, and zooplankton dynamics within the Caputo sense.The process of photosynthesis and the generation of oxygen in phytoplankton were studied by researchers 27 .They concentrated on the Caputo fractional derivative and compared the outcomes to the integerordered derivative.To understand how nutrients, hazardous phytoplankton, and zooplankton collaborate, a fractional-order mathematical model with a delay in time was developed in 28 .The major goal was to investigate how time latency and fractional order affected the ecology.Another study 29 used a fear function together with a Holling type II function to describe the various prey and predator organisms to explore the influence of memory on interactions between predators and prey in the context of global warming.The Caputo derivative and exponential decay function were used to examine how global warming affects both prey and predators.Bonyah 30 investigated a fractional-order model with two controls for nutrients, phytoplankton, and zooplankton.A broad formulation problem with state and adjoint equations, like a fractional optimal control issue.is provided in left fractional derivative terms.The forward-backward sweep approach, which is employed to address the fractional optimum control issue, was created using the Adams-type predictor-corrector technique.Kumar et al. 31 looked into an ecological model that includes the production of oxygen during photosynthesis, plankton apnea, and the effects of zooplankton predation on phytoplankton.This model is a fractional derivative and considers all three components of the food chain.Using a generalized Liouville-Caputo type fractional derivative, they have created a system of three non-integer order differential equations.They started by introducing the issue more conceptually and then employed a current fractional numerical method to demonstrate it experimentally.A fractional-order nonlinear mathematical model was put up in 32 to examine the behavior of climate change using the Caputo operator.The model explained the effects of declining water quality brought on by greenhouse gases on populations of aquatic animals.The model looked at equilibrium locations and talked about how equilibria can be asymptotically stable.The model's singular solution was established, and the numerical solution was discovered using a generalized predictor-corrector technique.As a result of the above debate, we analyze a fractional-order mathematical representation with time-dependent parameters to show the fast global warming's prospective impacts on marine ecosystems.The Caputo fractional derivative notion performs better than ordinary integer-order derivatives.This distinguishes our concept and process from the earlier model that has been provided up to this point.
A generalized version of the model and a summary of the description of the proposed model are provided in Section "Environment Management Model with Caputo fractional derivative".Furthermore, the theoretical background of the suggested fractional operator is explored.Section Qualitative analysis of the model" deals with the qualitative analysis of the proposed system.The numerical solutions to the suggested fractional-order model with power law kernel are provided in Section "Numerical scheme".In Sections "Numerical simulation" and "Conclusion, " the numerical simulations, results, and conclusions are addressed.

Environment management model with caputo fractional derivative
We take into consideration the ecological maintenance model with time-varying characteristics described in 5 .A model that explores how global warming would affect aquatic ecosystems as a result of the fast emitted greenhouse gases (GGs) caused by humankind.The four main components of the diverse system are the density of environmental GGs, G(t) , which are quickly released through different sources; rising temperature in the atmosphere, T(t) , which grows correspondingly with rising levels of environmental GGs along with is the cause of global warming; the density of planktonic population in marine ecosystems, P(t) , which is continually at risk by rising temperature and GGs concentration; and the fish population density in marine environments, F(t) , which is likewise declining in quantity because of increasing global warming, acidification, a lack of saturation oxygen, and a lack of planktonic species.

Model's assumptions
• The usual growth rates of G and T are g 1 and g 2 , respectively.However, in not having any of the negative effects of GGs and global warming, P and F expand at their usual rates of g 3 and g 4 , respectively.• In marine ecosystems, fish emit saturated Carbon-dioxide, which marginally raises the concentration of GGs.
The rise in GGs concentration caused by the fish population is shown in this case by ϕ 1 GF.• The photosynthesis carried out by phytoplankton in marine ecosystems lowers the levels of GGs in the atmosphere.Planktonic population GGs absorption is represented by ϕ 2 GP.• The density of GGs, defined by ϕ 3 T , rises as a result of natural disasters brought on by climate change, such as droughts and forest fires.• In accordance with the density of ambient GGs, the temperature of the atmosphere rises.A rise in the tem- perature of the atmosphere brought on by rising GGs is shown here in the form of 1 GT.• The photosynthesis of aquatic plankton is influenced by temperature, and this enables them to counteract increasing temperatures.Planktonic population absorption of ambient temperature is represented by the symbol 2 PT.• All living beings have an ultimate carrying capacity, therefore supposing a constant level, α(0 < α < 1) , dissolved Carbon-dioxide saturation is essential to their lifespan.• We choose Q 1 and Q 2 as the carrying capacities for the fish population and the planktonic population, respec- tively.
Q 2 are the appropriate decomposition rates.• While high concentrations impede development through increased plankton respiration and decreased oxy- gen dispersion density, they do enhance the overall density of marine plankton.This inhibits the growth of marine fisheries and causes a dissolved oxygen shortfall when paired with saturated Carbon-dioxide.As a result, a rise in planktonic population resulting from Carbon-dioxide absorption is represented by " ω 1 P α+G ′′ , while a fish population reduction caused by an excessive amount of dissolved Carbon-dioxide is shown by α+G .• Plankton and fish density are continuously declining as a result of acidification and rising temperatures harm- ing marine ecosystems.Planktonic population reduction caused by warming is shown by ω 2 PT , planktonic population reduction linked to acidity is shown by ω 4 GP , and the fish density demise in response to rising temperatures is shown by ξ 3 TF.• In marine ecosystems, fish populations serve as predators and planktonic populations as prey.Therefore, we infer that the decline in the planktonic population due to fish population predation is represented by ω 3 FP .And ξ 1 FP represents the growth in fish density caused by their consumption of planktonic organisms.
Here, we use the nonlinear fractional-order system to define the aforementioned description.
Vol:.( 1234567890) www.nature.com/scientificreports/Where C D β represents the Caputo derivative of order 0 < β ≤ 1 .The corresponding nonnegative initial condi- tions are such that Now we'll go through some recent and relevant calculus results.
Definition 2.1 33 The Caputo derivative of a differentiable function ζ(t) to order β ∈ (0, 1) with beginning point, t = 0 , is given by is an integrable function with 0 < β < 1 , the fractional integral is specified as follows 34 : Remark 2.1 A fixed point τ ⋆ is regarded as the equilibrium point of the Caputo system if and only if ζ(t, τ ⋆ ) = 0.

Well-posedness and Positively Invariant Region
Here, we examine the conditions necessary for a system to produce favorable results while taking into account realistic real-world scenarios.
Theorem 3.1 For every t ≥ 0 , the system (1)'s solutions are not negative.

Proof We have got
The choice of solution can't escape from the hyperplane if (G(0), T(0), P(0), F(0)) ∈ R 4 + .The vector field on each hyperplane surrounding the non-negative orthant directs into the domain R 4 + , making it a positively invariant set.

Equilibrium Points Analysis
If the system (1)'s left side is set to zero, the two different sorts of equilibrium points are obtained 5 are given here • The equilibrium point when there are no more fish is , www.nature.com/scientificreports/

Solutions' Existence and Uniqueness
Implementing the Banach fixed point theory and Schaefer's fixed point theorem, this section establishes the existence and distinctiveness of a solution to the system (1).Organize the following function: The Caputo fractional derivative model (1) of order β > 0 will therefore be subjected to the fractional integral according to its beginning conditions.The second kind of Volterra-integral equations that are produced as a consequence of the procedure provide the answer to the proposed model (1).
In order for (D, .) to be the Banach space and K 1 ([0, T]) to be the Banach space comprising all continuous functions established in where where where The fact that the condition (G, T, P, F) Ŵ(1−β) sin(βπ)T β βπ < 1 is clear from the data.Considering that the parameter J contains a fixed point in t ∈ [0, T] as it is a contraction mapping, the Banach contraction mapping concept is applied to demonstrate this.Now, we look at the existence of solutions for the system (1) using Schaefer's fixed point theorem. ( Vol:.(1234567890)

Scientific
where ρ is a random number between 0 and 1, then system (1) has at least one solution.

And
The above expressions approach zero when t a → t b on the right side of the inequality.According to the Arzela- Ascoli theorem, J is a continuous function.Now, we prove that is bounded for some 0 < µ < 1 by (1).For every t ∈ [0, T] , let (G, T, P, F) ∈ Q(J) , such that (G, T, P, F) = µJ(G, T, P, F) , yields

And
As we have established that Q(J) is bounded, the system (1)'s solution exists because J has a fixed point which is determined by Schaefer's fixed point theorem.

, T], D) if
Proof There exists ℑ = {δ, τ , φ, ς} > 0 such that is true for all t ∈ [0, T] according to definition (3.1), which designates P as a non-decreasing function of t.The continuous nature of the functions G , T , P , and F has been shown, and the Lipschitz condition is met when (G, T, P, F) > 0 .Theorem (3.3) provides a unique answer for the proposed system (1) When we integrate the inequalities in the definition (3.1), we acquire From equation (54) and Lemma (2.2), we have (49) Vol (56)

Numerical scheme
In the literature, it has been proposed that the power-law kernel-based Caputo derivative is suitable for mimicking power-law processes in practical problems.We use a numerical scheme based on a Newton polynomial to solve the system (1) numerically.

Numerical simulation
The model's numerical simulations have been carried out using the generalized two-step Lagrange polynomial for the power law kernel and the parametric values from 5 which are: δ 1 = 0.00095 , δ 2 = 0.0099 , δ 3 = 0.00025 , δ 4 = 0.0002 1000 , γ 1 = 0.0029 , γ 2 = 0.00099 , γ 3 = 1.0 , 1 = 0.00025 , 2 = 0.00565 , ω 1 = 0.00108 , ω 2 = 0.00001 , ω 3 = 0.0031 , ω 4 = 10.1 , ξ 1 = 175 1000 , ξ 2 = 190 1000 , ξ 3 = 61 1000 , α = 0.01 , Q 1 = 1000, 000 , and Q 2 = 10, 000 .We have used the value of �(t) = 0.01 .While the starting values of the continually changing species are G(0) = 0.04 , T(0) = 0.07 , P(0) = 17.5 , and F(0) = 7.8 .The simulations' main objectives are to confirm the analytical find- ings of this work and to clarify the dynamic behaviors of the organisms under consideration, particularly the plankton and fish populations in marine ecosystems under accelerated global warming.By dispersing energy, heat, and materials, marine ecosystems significantly contribute to restoring the equilibrium of the environment.Simulations of the proposed design model show that when analyzing internal behavior, the overall density of all the segments will fluctuate between 0 and 1. Analytical outcomes and to describe the model's application methodology examined the negative effects of the frequent, fast environmental concentration on marine ecosystems in Figs. 1, 2, 3 and 4 at fractional orders α = 1.0, 0.95, 0.85, 0.80 .Also, the impact of fractional order at different (78) K 4 (t q , G q , T q , P q , F q ) − 2K 4 (t q−1 , G q−1 , T q−1 , P q−1 , F q−1 ) Fractional-order derivations are more effective than traditional integer-order models in explaining physical processes.The present study examines the extent to which global warming will affect plankton and fish populations in marine bio-diversity, taking into account the impact of fractional memory.Additionally, the results at different fractal dimensions and bounded solitude in the domain shown in Figs.17, 18 and 19 are discussed, providing support for both theoretical and experimental observations.Additionally, this study forecasts the future of marine ecosystems, including fish and plankton populations, as well as rapid global warming through long-term numerical analysis of the dynamic behavior of dynamic organisms.

Conclusion
In this study, we proposed a fractional order environment management model to research how rapidly accumulating ambient greenhouse gases (GGs) are affecting marine ecosystems and how this is contributing to global warming.We investigated the positively invariant region and showed that the model has positive, limited solutions.This can shed light on the resilience of these ecosystems and help in devising strategies for their conservation.To examine the model's existence and uniqueness, we also applied methods from several fixed-point theorems.Our results show the model to be generalized Ulam-Hyers-Rassias stable.The mathematical model was then resolved using a numerical method based on Newton's polynomial interpolation.Elaborating on the specific numerical techniques employed and the accuracy of the results would provide a better understanding of the model's practical utility.Results using different fractional values show significant variations.Non-integer order greatly impacts the flexibility and behavior of the solution curves, as is visible from the graphs.It indicates that small perturbations in the model's parameters do not lead to drastic changes in behavior, which is relevant for decision-making.This research can be expanded to include more generalized applicable fractional operators and improve control strategies for marine fisheries resources.To optimally utilize marine fisheries resources,
[0, T] → D formed with Chebyshev norm, the functions (G, T, P, F) : [0, T] × D → D is considered to be continuous.The Lipschitz condition is met by the continuous functions G , T , P , and F if Therefore

Figure 17 .
Figure 17.Simulation of all compartments in feasible region with chaotic form at fractional order 1.

Figure 18 .
Figure 18.Simulation of all compartments in feasible region with chaotic form at fractional order 0.5.

Figure 19 .
Figure 19.Simulation of all compartments in feasible region with chaotic form and bounded solution at different fractional value.
β−1 dν consequence of fast global warming; however, fish populations decrease proportionately with a reduction in plankton diversity due to shortages of food.Thus, the growing amount of greenhouse gases encourages global warming, which drastically lowers the planktonic population through increased acidity and warming, and the associated outcomes are significantly reduced.A comparison of the results is drawn in Figs. 9, 10, 12, 13, 14, 15 and 16 by using power law, exponential law, and Mittag Leffler kernel at fractal dimension β = 0.9 and β = 0.8 respectively and solution bounded to the steady state point rapidly.Simulation of all compartments in feasible regions with chaotic form and bounded regions at different fractal fractional values is shown in Figs.17, 18 and  19.
Figure 2. Simulation of T(t) with Caputo fractional operator.valuesα=0.65,0.60, 0.55, 0.50 is shown in Figs.5, 6, 7 and 8)to observe the complete transmission by chang- ing values.The quantity of greenhouse gases in the natural environment is constantly rising, and this means that the atmospheric temperature is rising proportionately to the quick volume of GGs.This quick volume of GGs also leads to the introduction of acidification in ocean water, that eliminates both plankton diversity and fisheries resources in oceans.In addition, marine fisheries resources are in danger and reducing primarily as Figure 4. Simulation of F(t) with Caputo fractional operator.Vol.:(0123456789)Scientific Reports | (2023) 13:22441 | https://doi.org/10.1038/s41598-023-49806-7 Figure 7. Simulation of P(t) with Caputo fractional operator with changing fractional values .Figure 8. Simulation of F(t) with Caputo fractional operator with changing fractional values.a